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The NJL model has recently been extended with a method to simulate confinement. This leads in mean field 
approximation to a natural mechanism for the saturation of nuclear matter. We use the model to investigate the 
equation of state of asymmetric nuclear matter and then use it to compute the properties of neutron stars. 



1. Introduction 

The properties of isospin asymmetric matter 
at finite density are of particular interest in the 
study of nuclear physics and neutron stars. In 
particular, there is great interest in exploring 
whether stars with sufficiently high core density 
contain some form of quark matter |1I2I3| . In 
order to investigate such questions it is natural 
to start with a model for dense matter which is 
built from the quark level. A number of models of 
this kind have been constructed, for example the 
quark-meson coupling (QMC) model |4I5I6I7I8| 
and the chiral soliton model 0. 

We work within the framework of the flavor 
SU(2) Nambu-Jona-Lasinio (NJL) model, modi- 
fied to simulate quark confinement. As explained 
by Bentz and Thomas ^U], a description of nu- 
clear matter can be obtained, which naturally 
exhibits the property of saturation, thus allow- 
ing for stable nuclei at low densities. The nu- 
cleoli is constructed as a quark-diquark state by 
solving the Faddeev equation in the static ap- 
proximation Using the proper time cut-off 

on 



scheme to regularize the integrals |12I18| . effec- 
tively introduces the phenomena of confinement, 
while leading to the stabilization of symmetric nu- 
clear matter [TJ]|. The general form for the equa- 



tion of state is derived from the quark-level La- 
grangian using the path integral formalism |14| . 
This model has been developed further to incor- 
porate the phase transition to quark matter, in- 
cluding a color superconducting state ^3] How- 
ever, in the present work we will deal with central 
densities where such transitions have not yet oc- 
cured. 

In this work we consider isospin asymmetric 
matter (neutrons, protons and electrons) in (3- 
equilibrium. Electrons are introduced to the 
system as a Fermi gas, balancing the positive 
charge of the protons. We calculate neutron star 
masses, radii and profiles by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equations We 
find that this model can support stable stars up 
to a maximum mass of 2.19M where the central 
density is 0.97/to~ 3 . A typical mass neutron star 
(1.4M ©) has a central density of 0.4/m -3 , and a 
radius around 12.4 km. 



2. Model for the Equation Of State 

The Lagrangian for the flavor SU(2) NJL model 
has the general form 



C = WQf- m)i> + J2 Ga(^L Q V) 2 , 



(1) 
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Figure 1. Equation of state for asymmetric nu- 
clear matter for various values of the proton frac- 
tions (Yp). 



Figure 2. Proton fraction for charge-neutral nu- 
clear matter (neutrons, protons and electrons) in 
chemical equilibrium. 



where to is the current quark mass, ip is the flavor 
SU(2) quark field, T a are the matrices in Dirac, 
flavor and color space, characterizing the chi- 
rally symmetric 4-fermi interaction between the 
quarks, and G a are the coupling constants asso- 
ciated with these interactions. More explicitly, in 
matter of density p, we have 



C 



(M - to) 2 



4G W 



+ G^uj^ + G p p^ Pll + Ci , 



(2) 



where we define the effective quark mass (i.e., 
the constituent quark mass in-medium) as M = 
to — 2G n (p\ipip\p} and the vector and isovector 
fields in the nuclear medium are given by uj^ = 
{p\'tp r y ,J -ip\p) and p M = {p\$'f fJ, Tip\p). From this La- 
grangian we derive both the equation of state and 
the quark-diquark structure of the bound "nu- 
cleon" . The interaction Lagrangian £/ consists 
of a sum of qq channel interaction terms, which 
are isolated using Fierz transformations ^B] 



d = Cr+Cr + 



(3) 



In order to construct the nucleon state, we restrict 
ourselves to the scalar diquark channel (£/ s ) in 



the present calculation. This will soon be ex- 
tended to include vector di-quark correlations. 
The equation of state can be derived using the 
path integral formalism |14| . or equivalent La- 
grangian techniques |17j . The energy density is 
expressed as 



f — f 



£n + £uj + £ 



(4) 



where the vacuum term is given by 



= 12i 



d 4 k k 2 ~M 2 + ie (M-mf 



(2tt)4 
(M - to) 2 
4G T ' 



'k 2 ~M§ 



(5) 



The Fermi motion of the nucleons contributes 
d 3 k 



£ 



N 



N 



IN 



(2tt) : 



-9{k FN -k)^Jk 2 + M 2 N , (6) 



where Mm is the effective mass of the bound nu- 
cleon. (This is calculated by solving the Faddeev 
equations using the in-medium masses of both the 
constituent quark and scalar di-quark.) The sum 
here is over neutrons and protons (N = n,p), 
which give equal contributions in the case of sym- 
metric nuclear matter (of Fermi momentum kp N ). 
For nuclear matter at rest the omega meson con- 
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Figure 4. Neutron star masses as a function of 
central baryon density. 



Figure 3. Equations of state for neutron mat- 
ter (solid line), matter in (3 - equilibrium (dashed 
line) and symmetric nuclear matter (dotted line) . 



tribution is given by 



(7) 



where the baryon density ps, is the sum of the 
proton and neutron densities, p p and p n , respec- 
tively and G w is the vector coupling constant. 
The contribution to the energy density from the 
rho meson is given by, 



£ p = 9Gp(p p 



(8) 



where the coupling constant G p is calculated at 
the empirical saturation density, (pb, Eb / A) = 
(0.16/m -3 , — lhMeV) : so that the symmetry en- 
ergy coefficient is 32.5MeV. 

For a given density the constituent quark mass 
M in Eq. JSJ is determined by the condition 
d£/dM — 0, where the quark mass at zero den- 
sity is M = AOOMeV. The nucleon mass M N 
in Eq. 10 is a function of M, determined by 
the pole in the Faddeev equation [TJ]]. The en- 
ergy density and pressure are calculated implic- 
itly as a function of baryon density. As shown in 



Fig. ^ the presence of protons softens the equa- 
tion of state. This effect leads to more com- 
pact neutron stars. To find the equation of state 
for matter in /3-equilibrium, we introduce elec- 
trons, £ e — kp /47r 2 , and impose the conditions of 
charge neutrality (p p — p e ) and chemical equilib- 
rium (p e = p n — p p ). Fig. indicates the proton 
fraction for matter in /^-equilibrium as a function 
of baryon density. 



Table 1 

Parameters for the calculations of the properties 
of the nucleon and nuclear matter. 
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G a [GeV' 2 ] 


0.510^ 


G^GeV' 2 } 


0.37G W 


G p [GeV- 2 } 


O.C^G^ 


A uv [MeV] 


638.5 


A IR [MeV] 


200 


m[MeV] 


16.93 
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Figure 6. Pressure profiles for neutron stars of 
mass 1.4 MQ. 



Figure 5. Neutron star mass versus radius for 
pure neutron matter and for nuclear matter in 
/3-equilibrium. 



3. Results and Discussion 

Given the relationship between £ and P we can 
integrate the TOV equations, 

dP _ G{£{r) + P(r))(47rr 3 P(r) + M{r)) 

~d7 ~ r(r - 2GM(r)) ^ ' 

M = 4tt / £{r)r 2 dr (10) 
Jo 

The stars generated from matter in /?-equilibrium 
are similar to those composed of pure neutron 
matter. Fig. illustrates the profiles of typical 
mass stars for each equation of state. In the case 
of matter in /3-equilibrium the central density is 
a little less than 3 times normal nuclear matter 
density. This corresponds to a proton fraction of 
about 15% (see Fig.[2J), in the centre of the star. 

4. Concluding Remarks 

In conclusion, we have seen that matter in j3- 
cquilibrium produces more compact stars than 



pure neutron matter. The abundance of protons 
depends on the value of G p , which in our case 
is relatively low. While the differences are quite 
small, the effect is clear. We have found that this 
model predicts densities of about 0.4/to -3 in the 
core of typical mass neutron stars. Results for 
masses and radii are in basic agreement with as- 
trophysical data. 

This work was supported by the Australian 
Research Council and DOE contract DE-AC05- 
84ER40150, under which SURA operates Jeffer- 
son Lab. 
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